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Abstract. In popular cold dark matter cosmological scenarios, stars 
; may have first appeared in significant numbers around a redshift of 10 

. p. or so, as the gas within protogalactic halos with virial temperatures 

T v i r > 10 4 - 3 K (corresponding to masses comparable to those of present- 
day dwarf ellipticals) cooled rapidly due to atomic processes and frag- 
mented. It is this 'second generation' of subgalactic stellar systems, aided 
perhaps by an early population of accreting black holes in their nuclei, 
I \ which may have generated the ultraviolet radiation and mechanical en- 

^ ■ ergy that ended the cosmic "dark ages" and reheated and reionized most 

\ of the hydrogen in the universe by a redshift of 6. The detailed history 

of the universe during and soon after these crucial formative stages de- 
pends on the power-spectrum of density fluctuations on small scales and 
^ ■ on a complex network of poorly understood 'feedback' mechanisms, and 

(S( ! is one of the missing link in galaxy formation and evolution studies. The 

O ■ astrophysics of the epoch of first light is recorded in the thermal state, 

^ r*| ! ionization degree, and chemical composition of the intergalactic medium, 

CL| the main repository of baryons at high redshifts. 

o : 

in ■ 

^\ \ 1. Introduction 

> ■ 

At epochs corresponding to z ~ 1000 the intergalactic medium (IGM) is ex- 
pected to recombine and remain neutral until sources of radiation and heat 
develop that are capable of reionizing it. The detection of transmitted flux 
shortward of the Lja wavelength in the spectra of sources at z > 5 implies that 
the hydrogen component of this IGM was ionized at even higher redshifts. The 
increasing thickening of the Lya forest recently measured in the spectra of SDSS 
z ~ 6 quasars (Becker et al. 2001; Djorgovski et al. 2001) may be the signature 
of the trailing edge of the cosmic reionization epoch. It is clear that substantial 
sources of ultraviolet photons and mechanical energy, like young star-forming 
galaxies, were already present back then. The reionization of intergalactic hydro- 
gen at z > 6 is unlikely to have been accomplished by quasi-stellar sources: the 
observed dearth of luminous optical and radio-selected QSOs at z > 3 (Shaver 
et al. 1996; Fan et al. 2001), together with the detection of substantial Lyman- 
continuum flux in a composite spectrum of Lyman-break galaxies at (z) = 3.4 
(Steidel et al. 2001), both may support the idea that massive stars in galactic 
and subgalactic systems - rather than quasars - reionized the hydrogen com- 
ponent of the IGM when the universe was less than 5% of its current age, and 
dominate the 1 ryd metagalactic flux at all redshifts greater than 3. An episode 
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of pregalactic star formation may also provide a possible explanation for the 
widespread existence of heavy elements (like carbon, oxygen, and silicon) in the 
IGM. There is mounting evidence that the double reionization of helium may 
have occurred later, at a redshift of 3 or so (see Kriss et al. 2001, and references 
therein): this is likely due to the integrated radiation emitted above 4 ryd by 
QSOs. 

Establishing what ended the dark ages and when is important for deter- 
mining the impact of cosmological reionization and reheating on several key 
cosmological issues, from the role reionization plays in allowing protogalactic 
objects to cool and make stars, to determining the thermal state of baryons at 
high redshifts and the small-scale structure in the temperature fluctuations of 
the cosmic microwave background. Conversely, probing the reionization epoch 
may provide a means for constraining competing models for the formation of cos- 
mic structures: for example, popular modifications of the CDM paradigm that 
attempt to improve over CDM by suppressing the primordial power-spectrum 
on small scales, like warm dark matter (WDM), are known to reduce the num- 
ber of collapsed halos at high redshifts and make it more difficult to reionize 
the universe (Barkana et al. 2001). In this talk I will summarize some recent 
theoretical developments in understanding the astrophysics of the epoch of first 
light and the impact that some of the earliest generations of stars, galaxies, and 
black holes in the universe may have had on the IGM. 



2. Reionization by massive stars 

In a cold dark matter (CDM) universe, structure formation is a hierarchical 
process in which non linear, massive structures grow via the merger of smaller 
initial units. Large numbers of low-mass galaxy halos are expected to form at 
early times in these popular theories, leading to an era of reionization, reheating, 
and chemical enrichment. Most models predict that intergalactic hydrogen was 
reionized by an early generation of stars or accreting black holes at z = 7 — 
15. One should note, however, that while numerical N-body+hydrodynamical 
simulations have convincingly shown that the IGM does fragment into structures 
at early times in CDM cosmogonies (e.g. Cen et al. 1994; Zhang et al. 1995; 
Hernquist et al. 1996), the same simulations are much less able to predict the 
efficiency with which the first gravitationally collapsed objects lit up the universe 
at the end of the dark age. 

The scenario that has received the most theoretical studies is one where 
hydrogen is photoionized by the UV radiation emitted either by quasars or by 
stars with masses > 10 M (e.g. Shapiro & Giroux 1987; Haiman & Loeb 1998; 
Madau et al. 1999; Chiu & Ostriker 2000; Ciardi et al. 2000), rather than ionized 
by collisions with electrons heated up by, e.g. supernova-driven winds from 
early pregalactic objects. In the former case a high degree of ionization requires 
about 13.6 x (l+t/t rec ) eV per hydrogen atom, where i rec is the volume-averaged 
hydrogen recombination timescale, t/t rec being much greater than unity already 
at z ~ 10 according to numerical simulations. Collisional ionization to a neutral 
fraction of only a few parts in 10 5 requires a comparable energy input, i. e. an 
IGM temperature close to 10 5 K or about 25 eV per atom. 
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Massive stars will deposit both radiative and mechanical energy into the 
interstellar medium of protogalaxies. A complex network of 'feedback' mecha- 
nisms is likely at work in these systems, as the gas in shallow potential is more 
easily blown away (Dekel & Silk 1986; Tegmark et al. 1993; Mac Low & Ferrara 
1999; Mori et al. 2002) thereby quenching star formation. Furthermore, as the 
blastwaves produced by supernova explosions reheat the surrounding intergalac- 
tic gas and enrich it with newly formed heavy elements (see below), they can 
inhibit the formation of surrounding low-mass galaxies due to 'baryonic strip- 
ping' (Scannapieco et al. 2002). It is therefore difficult to establish whether an 
early input of mechanical energy will actually play a major role in determining 
the thermal and ionization state of the IGM on large scales. What can be easily 
shown is that, during the evolution of a a 'typical' stellar population, more en- 
ergy is lost in ultraviolet radiation than in mechanical form. This is because in 
nuclear burning from zero to solar metallicity (Z@ = 0.02), the energy radiated 
per baryon is 0.02 x 0.007 x mnc 2 , with about one third of it going into H-ionizing 
photons. The same massive stars that dominate the UV light also explode as 
supernovae (SNe), returning most of the metals to the interstellar medium and 
injecting about 10 51 ergs per event in kinetic energy. For a Salpeter initial mass 
function (IMF), one has about one SN every 150 M of baryons that forms 
stars. The mass fraction in mechanical energy is then approximately 4 x 10 -6 , 
ten times lower than the fraction released in photons above 1 ryd. 

The relative importance of photoionization versus shock ionization will de- 
pend, however, on the efficiency with which radiation and mechanical energy 
actually escape into the IGM. Consider, for example, the case of an early gener- 
ation of subgalactic systems collapsing at redshift 9 from 2-cr fluctuations. At 
these epochs their dark matter halos would have virial radii r v i r = 0.75/i -1 Kpc 
and circular velocities V c (r v i T ) = 25 kms -1 , corresponding in top-hat spherical 
collapse to a virial temperature T V [ T = 0.5fim p V^ /k « 10 4 3 K and halo mass 
M = G.IV^/GH « 10 8 /i -1 M .[] Halos in this mass range are characterized 
by very short dynamical timescales (and even shorter gas cooling times due to 
atomic hydrogen) and may therefore form stars in a rapid but intense burst 
before SN 'feedback' quenches further star formation. For a star formation ef- 
ficiency of / = 0.1, fifr/i 2 = 0.020, h = 0.5, and a Salpeter IMF, the explosive 
output of 10, 000 SNe will inject an energy Eq « 10 55 ergs. This is roughly a 
hundred times higher than the gas binding energy: a significant fraction of the 
halo gas will then be lifted out of the potential well ('blow-away') and shock 
the intergalactic medium (Madau et al. 2001). If the explosion occurs at cosmic 
time t = 4 x 10 8 yr (corresponding in the adopted cosmology to z = 9), at time 
At = OAt after the event it is a good approximation to treat the cosmological 
blast wave as adiabatic, with proper radius given by the standard Sedov-Taylor 
self-similar solution, 

fl ^(i^y'> A(2/SKl7kpc (1) 



lr rhis assumes an Einstein-de Sitter universe with Hubble constant Ho = 100 h kms 1 Mpc 1 . 

2 Here Qb is the baryon density parameter, and fQb is the fraction of halo mass converted into 
stars. 
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Here 77 « 0.3 is the fraction of the available SN energy that is converted into 
kinetic energy of the blown-away material (Mori et al. 2002). At this instant 
the shock velocity relative to the Hubble flow is 

v s « 2R s /5At « 40 kms~\ (2) 

lower than the escape velocity from the halo center. The gas temperature just 
behind the shock front is T s = 3fim p v 2 s /16k > 10 4 3 K, enough to ionize all in- 
coming intergalactic hydrogen. At these redshifts, it is the onset of Compton 
cooling off cosmic microwave background photons that ends the adiabatic stage 
of blast wave propagation; the shell then enters the snowplough phase and is 
finally confined by the IGM pressure. According to the Press-Schechter formal- 
ism, the comoving abundance of collapsed dark halos with mass M = 10 8 h~ x M 
at z — 9 is dn/dlnM ~ 80 /i 3 Mpc -3 , corresponding to a mean proper distance 
between neighboring halos of ~ 15 kpc. With the assumed star formation 
efficiency, only a small fraction, about 4 percent, of the total stellar mass in- 
ferred today (Fukugita et al. 1998) would actually form at these early epochs. 
Still, our simple analysis shows that the blast waves from such a population of 
pregalactic objects could drive vast portions of the IGM to a significantly higher 
adiabat, Tigm ~ 10 5 K, than expected from photoionization, so as to 'choke off' 
the collapse of further M < 10 9 M systems by raising the cosmological 
Jeans mass. In this sense the process may be self-regulating. 

The thermal history of expanding intergalactic primordial gas at the mean 
density is plotted in Figure 1 as a function of redshift for a number of illustrative 
cases. The code we have used includes the relevant cooling and heating processes 
and follows the non-equilibrium evolution of hydrogen and helium ionic species 
in a cosmological context. 




redshift 

Fi gure 1. Thermal history of intergalactic gas at the mean density in an 
Einstein-de Sitter universe with Q5/1 2 = 0.02 and h = 0.5. Short dash-dotted 
line: temperature evolution when the only heating source is a constant ultravi- 
olet (CUV) background of intensity 10 -22 ergs cm -2 s _1 Hz -1 sr _1 at 1 Ryd and 
power-law spectrum with energy slope a = 1. Long dash-dotted line: same for 
the time-dependent quasar ionizing background as computed by Haardt & Madau 
(1996; HM). Short-dashed line: heating due to a CUV background but with an 
initial temperature of 4 x 10 4 K at z = 9 as expected from an early era of pregalactic 
outflows. Long-dashed line: same but for a HM background. Solid line: heating 
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due to a CUV background but with an initial temperature of 10 5 K at z = 9. Dotted 
line: same but for a HM background. 

The gas is allowed to interact with the CMB through Compton cooling and ei- 
ther with a time-dependent QSO ionizing background as computed by Haardt 
& Madau (1996) or with a time-independent metagalactic flux of intensity 
10 -22 ergs cm -2 s _1 Hz -1 sr _1 at 1 Ryd (and power-law spectrum with energy 
slope a = 1). The temperature of the medium at z = 9 - where we start our 
integration - has been either computed self-consistently from photoheating or 
fixed to be in the range 10 4 6 — 10 5 K expected from SN-driven bubbles with 
significant filling factors (see below). The various curves show that the tem- 
perature of the IGM at z = 3 — 4 will retain little memory of an early era of 
pregalactic outflows and preheating, and be consistent with that expected from 
photoionization. 

3. Feedback and pregalactic enrichment 

Understanding the origin of the chemical elements, following the increase in their 
abundances with cosmic time, and uncovering the processes responsible for dis- 
tributing the products of stellar nucleosynthesis over very large distances are all 
key aspects of the evolution of gaseous matter in the universe. One of the major 
discoveries with Keck concerning the IGM has been the identification of metal 
absorption lines associated with many of the Lya forest systems. The detection 
of measurable C iv and Si iv absorption lines in clouds with H I column densi- 
ties as low as 10 14 5 cm -2 implies a minimum universal metallicity relative to 
solar in the range [—3.2] to [—2.5] at z ~ 3 (Songaila 1997). There is no indica- 
tion in the data of a turnover in the C iv column density distribution down to 
N C w £ 10 11 - 7 cm" 2 (JVri < 10 14 - 2 cm" 2 , Ellison et al. 2000). 

From a theoretical perspective, it is unclear whether the existence of heavy 
elements in the IGM at z = 3 — 3.5 points to an early (z > 6) enrichment epoch 
by low-mass subgalactic systems (Madau et al. 2001), or is rather due to late 
pollution by the population of star-forming galaxies known to be already in 
place at z = 3. The Press-Schechter theory for the evolving mass function of 
dark matter halos predicts a power-law dependence, dN/dlnm oc m( neff-3 )/ 6 , 
where n e ff is the effective slope of the CDM power spectrum, n e ff ~ —2.5 on sub- 
galactic scales. As hot, metal-enriched gas from SN-driven winds escapes its 
host halo, shocks the IGM, and eventually forms a blast wave, it sweeps a region 
of intergalactic space the volume of which increases with the 3/5 power of the 
injected energy Eq (in the adiabatic Sedov-Taylor phase). The total fractional 
volume or porosity, Q, filled by these 'metal bubbles' per unit explosive energy 
density E dN/d\nm is then Q oc dN/dlnm oc (dN/dlnm) 2 / 5 oc m" 11 / 30 . 
Within this simple scenario it is the star-forming objects with the smallest 
masses which will arguably be the most efficient pollutant of the IGM on large 
scales. Metal-enriched material from SN ejecta may be far more easily acceler- 
ated to escape velocities in the shallow potential wells of subgalactic systems at 
z > 6. Late enrichment may also encounter problems in explaining the kinematic 
quiescence of C iv lines; the observed small scale properties of the IGM at z ~ 3 
appear consistent with C iv absorbers being the result of ancient pregalactic 
outflows (Rauch et al. 2001). According to numerical hydrodynamics Simula- 
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tions of structure formation in the IGM, the metals associated with log TVhi < 14 
filaments fill a fraction > 5% of intergalactic space, and are therefore far away 
from the high overdensity peaks where galaxies form, gas cools, and star forma- 
tion takes place. Their chemical enrichment may then reflect more uniform (i.e. 
'early') rather than in-situ (i.e. 'late') metal pollution. The case for pregalactic 
enrichment may have been recently strenghtened by the observation of an invari- 
ant C IV column density distribution throughout the redshift range 2 < z < 5 
(Songaila 2001). 
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Fi gure 2. Snapshots of the logarithmic number density of the gas at five differ- 
ent elapsed times for our Case 1 simulation run. The three panels in each row show 
the spatial density distribution in the X — Y plane on the nested grids. The three 
columns in each figure depict the time evolution from about 6 Myr to up to 180 
Myr. Along a given row, the leftmost panel refers to grid L5 (linear size 48 kpc), 
the central one to grid L3 (linear size 12 kpc), and the rightmost panel refers to the 
grid LI (linear size 3 kpc). The density range is —5 < log (n/cm -3 ) < 1. The halo 
gas is assumed to be initially in hydrostatic equilibrium and non self-gravitating. 
(From Mori et al. 2002.) 
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To simulate the process of blow-away we have developed a 3D Eulerian code that 
solves the hydrodynamic equations for a perfect fluid in Cartesian geometry. To 
deal with very different length scales in our simulation we have adopted a 'nested 
grid method' with six levels of fixed Cartesian grids. The grids are connected 
by the transfer of conserved variables, and are centered within each other, with 
the finest covering the whole galaxy halo. Since the cell number is the same 
(128 x 128 x 128) for every level Ln (n = 1, 2, 6), the minimum resolved scale 
is about 22 pc and the size of the coarsest grid is 96 kpc. Thus, the scheme has 
a wide dynamic range in the space dimension. 

The results of a numerical simulation (run on massive parallel supercomput- 
ers at the Center for Computational Physics, Tsukuba University) are depicted 
in Figure 2. In this run about 10,000 SNe explode in a high redshift protohalo 
(corresponding to a star formation efficiency of about 10% for a Salpeter IMF). 
Our algorithm for simulating SN feedback improves upon previous treatments 
in several ways. OB associations are distributed as a function of gas density 
according to a Schmidt-type law (oc p a ) using a Monte-Carlo procedure. After 
a main sequence lifetime, all stars more massive than 8 M explode instanta- 
neously injecting an energy of 10 51 ergs, and their outer layers are blown out 
leaving a compact remnant of 1.4 M . Therefore SNe inject energy (assumed in 
pure thermal form) and mass into the interstellar medium: these are supplied 
to a sphere corresponding to the radius of a SN remnant in a uniform ambient 
medium of density n and in the adiabatic Sedov-Taylor phase; the expansion ve- 
locity is also self-consistently calculated from such solution. The gas then starts 
cooling immediately according to the adopted gas cooling function. While the 
time-step At is controlled by the Courant condition, if there are more than 
two SN events in a OB association we decrease At until only one explosion per 
association occurs during the time-step. 

Snapshots of the gas density distribution in a planar slice of the nested grids 
are shown for an extended stellar distribution case (a = 1). After a few Myr 
from the beginning of the simulation the most massive stars explode as SNe and 
produce expanding hot bubbles surrounded by a cooling dense (n ~ 1 cm -3 ) 
shell. As the evolution continues, a coherent and increasingly spherical shell 
expanding into the IGM is eventually formed. The shell contains a large fraction 
of the halo gas that has been swept-up during the evolution. The final bottom 
row of the simulation figure shows the final stages of the evolution: the shell is 
now nearly spherically symmetric, its interior being filled with warm (T < 10 6 K) 
gas at a very low density n < 10 -4 cm -3 , i.e. slightly below the mean value 
for the IGM. At the end of the simulations, the shell is still sweeping out IGM 
material; its radius and velocity are 21 kpc and 26 km sec -1 at t — 250 Myr. 
Using momentum conservation one can estimate the final radius of the shell to 
be close to 25 kpc. This is about 17 times the virial radius of the halo! 

It is clear then that SN-driven pregalactic outflows may be an efficient 
mechanism for spreading metals around. The collective explosive output of 
about ten thousands SNe per M > 10 8 h~ x M halo at these early epochs could 
then pollute the entire intergalactic space to a mean metallicity (Z) = Qz/^b 
0.003 (comparable to the levels observed in the Lya forest at z ~ 3) without 
much perturbing the IGM hydrodynamically, i.e. producing large variations of 
the baryons relative to the dark matter. The significance of these results goes 
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well beyond early metal enrichment. This is because, since the cooling time of 
collisionally ionized high density gas in small halos at high redshifts is much 
shorter than the then Hubble time, virtually all baryons are predicted to sink 
to the centers of these halos in the absence of any countervailing effect (White 
&; Rees 1978). Efficient feedback is then necessary in hierarchical clustering 
scenarios to avoid this 'cooling catastrophe', i.e. to prevent too many baryons 
from turning into stars as soon as the first levels of the hierarchy collapse. The 
required reduction of the stellar birthrate in halos with low circular velocities 
may naturally result from the heating and expulsion of material due to OB 
stellar winds and repeated SN explosions from the first burst of star formation. 



4. Probing the epoch of first light at 21— cm 

Prior to the epoch of full reionization, the intergalactic medium and gravita- 
tionally collapsed systems may be detectable in emission or absorption against 
the CMB at the frequency corresponding to the redshifted 21-cm line (associ- 
ated with the spin-flip transition from the triplet to the singlet state of neutral 
hydrogen.) In general, 21-cm spectral features are expected to display angular 
structure as well as structure in redshift space due to inhomogeneities in the gas 
density field, hydrogen ionized fraction, and spin temperature. Several different 
signatures have been investigated in the recent literature: (a) the fluctuations 
in the redshifted 21-cm emission induced by the gas density inhomogeneities 
that develop at early times in CDM-dominated cosmologies (Madau et al. 1997; 
Tozzi et al. 2000) and by virialized "minihalos" with T vir < 10 4 K (Iliev 2002); 
(b) the sharp absorption feature in the radio sky due to the rapid rise of the 
hya continuum background that marks the birth of the first UV sources in the 
universe (Shaver et al. 1999); (c) the 21-cm narrow lines generated in absorption 
against very high-redshift radio sources by the neutral IGM (Carilli et al. 2002) 
and by intervening minihalos and protogalactic disks (Furlanetto & Loeb 2002). 

A quick summary of the physics of 21-cm radiation will illustrate the ba- 
sic ideas and unresolved issues behind these studies. The emission or absorp- 
tion of 21-cm photons from a neutral IGM is governed by the hydrogen spin 
temperature T$ defined by ni/no = 3 exp(— T^/T^), where no and n\ are the 
singlet and triplet n — 1 hyperfine levels and fcgT* = 5.9 x 10 -6 eV is the en- 
ergy of the 21-cm transition. In the presence of only the CMB radiation with 
Tqmb = 2.73 (1 + z) K, the spin states will reach thermal equilibrium with the 
CMB onatimescale of T*/(Tcmb-4io) ~ 3xl0 5 (l+^) _1 yr (^io = 2.9xl0~ 15 s- 1 
is the spontaneous decay rate of the hyperfine transition of atomic hydrogen), 
and intergalactic H I will produce neither an absorption nor an emission signa- 
ture. A mechanism is required that decouples T# and Tcmb, e.g. by coupling 
the spin temperature instead to the kinetic temperature Tk of the gas itself. 
Two mechanisms are available, collisions between hydrogen atoms (Purcell & 
Field 1956) and scattering by hya photons (Field 1958). The collision-induced 
coupling between the spin and kinetic temperatures is dominated by the spin- 
exchange process between the colliding hydrogen atoms. The rate, however, is 
too small for realistic IGM densities at the redshifts of interest, although col- 
lisions may be important in dense regions with Sp/p > 30[(1 + z)/10] -2 , like 
virialized minihalos. 
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In the low density IGM, the dominant mechanism is the scattering of con- 
tinuum UV photons redshifted by the Hubble expansion into local Lya photons. 
The many scatterings mix the hyperfine levels of neutral hydrogen in its ground 
state via intermediate transitions to the 2p state, the Wouthuysen-Field process. 
As the neutral IGM is highly opaque to resonant scattering, the shape of the 
continuum radiation spectrum near Lya will follow a Boltzmann distribution 
with a temperature given by the kinetic temperature of the IGM (Field 1959). 
In this case the spin temperature of neutral hydrogen is a weighted mean be- 
tween the matter and CMB temperatures. There exists then a critical value of 
the background flux of Lya photons which, if greatly exceeded, would drive the 
spin temperature away from Tcmb- 

While the microphysics is well understood, our understanding of the as- 
trophysics of 21-cm tomography is still poor. In Tozzi et al. (2000) we used 
N-body cosmological simulations and, assuming a fully neutral medium with 
Ts ^> Tcmb, showed that prior to reionization the same network of sheets and 
filaments (the 'cosmic web') that gives rise to the Lya forest at z ~ 3 should 
lead to fluctuations in the 21-cm brightness temperature at higher redshifts 
(Figure 3). At 150 MHz (z = 8.5), for observations with a bandwidth of 1 
MHz, the root mean square fluctuations should be ~ 10 mK at 1', decreas- 
ing with scale. Because of the smoothness of the CMB sky, fluctuations in 
the 21-cm radiation will dominate the CMB fluctuations by about 2 orders of 
magnitude on arcmin scales. The search at 21-cm for the epoch of first light 
has become one of the main science drivers of the LOw Frequency ARray (LO- 
FAR; see |htt p : / / www . astron . nl /lof ar / science / ) . While remaining an extremely 
challenging project due to foreground contamination from extragalactic radio 
sources (Di Matteo et al. 2002), the detection and imaging of these small-scale 
structures with LOFAR is a tantalizing possibility within range of the thermal 
noise of the array. 
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Figure 3. Left: Radio map of redshifted 21-cm emission against the CMB in a 
'tilted' CDM (tCDM) cosmology at z — 8.5. Here a collisionless iV-body simulation 
with 64 3 particles has been performed with Hydra (Couchman et al. 1995). The 
simulation box size is 20h~ 1 comoving Mpc, corresponding to 17 (11) arcmin in 
tCDM (OCDM). The baryons are assumed to trace the dark matter distribution 
without any biasing, the spin temperature to be much greater than the temperature 
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of the CMB everywhere, and the gas to be fully neutral. The point spread function 
of the synthesized telescope beam is a spherical top-hat with a width of 2 arcmin. 
The frequency window is 1 MHz around a central frequency of 150 MHz. The 
contour levels outline regions with signal greater than 4 /iJy per beam. Right: Same 
for a open cosmology (OCDM). Since the growth of density fluctuations ceases early 
on in an open universe (and the power spectrum is normalized to the abundance of 
clusters today), the signal at a given angular size is much larger in OCDM than in 
tCDM at these early epochs. (From Tozzi et al. 2000.) 



On the theoretical side, there are several effects that need to be examined fur- 
ther. As mentioned above, it is the presence of a sufficient flux of Lja photons 
which renders the neutral IGM 'visible'. Without heating sources, the adia- 
batic expansion of the universe will lower the kinetic temperature of the gas 
well below that of the CMB, and the IGM will be detectable through its ab- 
sorption. If there are sources of radiation that preheat the IGM, it may be 
possible to detect it in emission instead. The energetic demand for heating the 
IGM above the CMB temperature is meager, only ~ 0.004 eV per particle at 
z ~ 10. Consequently, even relatively inefficient heating mechanisms may be 
important warming sources well before the universe was actually reionized. Per- 
haps more importantly, prior to full reionization the IGM will be a mixture of 
neutral, partially ionized, and fully ionized structures: low-density regions will 
be fully ionized first, followed by regions with higher and higher densities. Radio 
maps at 21 (1 + z) cm will show a patchwork of emission/absorption signals from 
H I zones modulated by H n regions where no signal is detectable against the 
CMB. It is the early generation of stars likely responsible for reionization which 
will also generate a background radiation field of photons with energies between 
10.2-13.6 eV to which the IGM is transparent. As each of these photons gets 
redshifted, it ultimately reaches the Lja transition energy of 10.2 eV, scatters 
resonantly off neutral hydrogen, and mixes the hyperfine levels. The observ- 
ability of the pre-reionization IGM depends in this case on the UV spectrum 
of the firsts stars, i.e. on the number of Lja continuum photons emitted per 
H-ionizing photon. 




Figure 4. Left: Continuum Lya photon emission rate N a (s 1 ) versus age (in 
years) for a stellar population with a Salpeter IMF (dashed lines) and a Scalo IMF 
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(solid lines). Upper curves assume an instantaneous burst of star formation, lower 
curves a constant star formation rate. The rate of production of Lya photons has 
been normalized to the rate of emission of hydrogen-ionizing photons, iVion- The 
population synthesis values have been computed for low metallicities (Z = 0.02 Zq), 
and are based on an update of Bruzual & Chariot's (1993) libraries. Right: Same 
after integration over the age of the stellar population. 



The integrated UV spectrum of a stellar population is characterized by a strong 
break at the Lyman edge whose size depends on age, initial mass function, and 
star formation history. In Figure 4 we have therefore normalized the Lyo: rate to 
the rate of emission of hydrogen-ionizing photons; after integrating over the stel- 
lar age, r, a starburst is characterized by a ratio N a /N[ on « 22 (r/0.1 Gyr) 0-22 
(Scalo IMF), about three times higher than the constant SFR case (this is be- 
cause of the shorter lifetime of the massive stars that produce 1 ryd radiation). 
It is possible that the startburst mode may be more relevant for the galaxies 
responsible for reionization, as stellar winds and supernovae can easily expell 
the gas out of the shallow potential wells of low mass halos, after the first burst 
of star formation occurs. In this case the escape fraction might be low until 
the gas is expelled, and as Figure 4 indicates this would contribute to an even 
stronger Lya continuum emission per ionizing flux. 
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